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THRESHOLD  COMPARISON  OF  PHASE-LOCK,  FREQUENCY-LOCK  AND 
MAXIMUM- LIKELfflOOD  TYPES  OF  FM  DISCRIMINATORS  * 

J.  J.  Spllker,  Jr. 

Communications  and  Controls  Research 
Lockheed  Missiles  and  Space  Co. 

Sunnyvale,  California 


Summary  -  In  making  a  comparison  of  demodulation 
techniques  for  FM  signals,  one  of  the  most  important 
criteria  for  evaluation  Is  the  threshold  input  signal - 
to-nolse  ratio  (SNR).  It  Is  at  this  value  of  input  SNR 
that  the  output  SNR  begins  to  decrease,  substantially 
more  rapidly  than  a  simple  proportionality  to  the  input 
SNR. 

The  objective  of  this  paper  la  to  present  a  consistent 
derivation  of  the  threshold  SNR  for  several  major  types 
of  FM  discriminators,  namely,  the  phase-lock  dis¬ 
criminator,  the  frequency-lock  discriminator  (FM  dis¬ 
criminator  with  negative  feedback),  and  the  maximum- 
likelihood  (a  posteriori,  most  probable)  computer.  This 
last  discriminator  Is  of  Interest  mainly  because  Its  op¬ 
eration  provides  a  bound  on  the  performance  of  other 
tj’pes. 

The  operation  of  the  discriminators  in  the  presence  of 
gausslan  Interfering  noise  Is  described,  and  the  causes 
of  threshold  are  investigated.  The  probability  of  losing 
the  "locked  on"  state  of  operation  is  computed.  Curves 
of  the  threshold  Input  SNR  are  plotted  vs.  the  ratio  of 
signal -to-modulatlon  bandwldths,  and  these  curves  are 
compared  for  the  various  types  of  discriminators.  Ex¬ 
perimental  confirmations  of  some  of  the  theoretical 
results  are  given. 

Introduction 

In  nniklng  a  comparison  of  FM-dcmodid-ition  lechniques, 
one  of  the  most  important  criteria  for  evaluation  is  the 
threshold  Input  SNR.  It  is  at  this  value  of  input  SNR 
that  the  output  SNR  Ijeglnsto  decrease,  substantiallv 
more  rapidly  than  a  simple  proportlonalltv  to  the  input 
SNR.  Wien  a  discriminator  is  operating  well  .above  Its 
threshold  and  utilizes  optimum  filtering,  iti  an  beshmvn 
that  the  output  SNR  is  linearly  related  to  Ihc  input  SNR 
by  an  expression  that  applies  to  almost  any  typo  of  dis¬ 
criminator. 

The  main  contribution  of  this  paper  Is  the  deriv.ation 
and  comparison  of  the  lhre..Jhold  Input  SNR  for  several 
major  types  of  KM  disorimln.ators,  namely,  the  phase- 
lock  discriminator,  the  frcciuencv-lock  di.scriminalor 
(FM  discriminator  with  negative  feedback),  and  the 
m.Tximum -likelihood  (a  posteriori,  mo.st  probable) com¬ 
puter.  The  operation  .and  threshold  behavior  of  these 
discriminators  is  analyzed  in  some  detail,  and  tholoop 
filters  are  optimized  in  order  to  minimize  the  v.alue  of 
the  threshold  input  SNR. 


*This  work  was  supported  in  part  in  ARDC.  Rome  ADC, 
Criffiss  AFn.  N.  V. 


Of  the  techniques  discussed  In  this  paper,  an  early 
version  of  the  frequency-lock  technique  was  the  first 
to  be  presented  In  the  literature.  In  19.^9, J.  G.  Chaffee*' 
described  the  use  of  feedback  in  an  FM -demodulation 
operation.  This  discriminator  did  not  require  a  lim¬ 
iter,  but  served  both  as  an  alternative  to  a  limiter  and 
as  a  method  of  reducing  distortion.  1, 2, 3, 4  j^ore 
recently,  interest  in  the  frequency-lock  technique  has 
been  revived  as  a  means  for  reducing  the  threshold 
Input  SNR,  ^  However,  to  the  author's  knowledge,  no 
analysis  of  the  threshold  improvement  Is  available  In 
the  ]  Herat  me. 

The  frequency-lock  discriminator  an.alyz.ed  In  this 
paper  differs  from  the  one  used  by  Chaffee  in  that  an 
ideal  frequency  discriminator  Is  employed  inside  the 
feedback  loo|),  whereas  Chaffee's  discriminator  didnot 
use  a  iiniltcr  and  w,ns  amplitude  sensitive.  (,4n  "ideal" 
frequency  discriminator  Is  defined  ns  a  ilevlce  where 
Vvdtage  output  is  proportional  to  the  InsHmtaneous  fre¬ 
quency  of  Its  input.  It  can  be  accurately  approximated 
l)V  a  wide-h.".rnl  limiter  FM  di.scriininator ,  or  z.ero- 
crossing  eomiting  device.)  This  paper  analyzes  the 
performance  and  causes  of  threshold  in  the  frequency- 
loek  ili.seriminalor  ami  optimizes  the  loop  fillers  so  n.s 
to  minimize  the  Ihreshold  input  SNR. 

rhe  phnsc-iock  FM  discriminator  and  other  phase- 
locked  circuits  h.ive  received  considerably  more  atten¬ 
tion  in  the  lller.nture  than  has  the  frctiucncy-lock  dis¬ 
criminator.  The  first  p.nprrs  analyzing  this  technique 
.seem  to  have  appeared,  about  19.')1,  <>»  < .  8, 9  This  set 
of  papers  dis'-usses  phase -locked  circuits  as  applied  to 
the  problem  of  recovering  the  reference  phase  of  a 
color  TV  signal.  One  of  the  more  Important  appllca- 
iions  of  phase-ioeked  rlrcuifs  is  in  tracking  .systems, 

.and  numerous  papers  have  appeared  in  the  literature 
dealing  with  the  design  of  phase-locked  circuits  optim¬ 
ized  to  tr.ack  .appropriate  transients  of  phase.  10, 11, 12 
Other  [xapers  have  presented  new  design  .and  analysis 
techniques.  13,  14 

This  treatment  of  the  phase-lock  discriminator  differs 
from  previous  work  in  that  ii  i.s  aimed  specifically  at 
the  problem  of  demodulating  an  FM  .signal  which  i.s 
modul.ated  by  analog  intelligence  having  a  low-pass 
spectrum.  The  purpose  of  the  analysis  is  to  determine 
the  threshold  behavior  of  a  phase-lock  discriminator 
which  ha.s  lieen  optimized  (in  a  restricted  .sense)  for 
this  type  rif  signal . 

rhe  maximuni  i i'Keitiieed  k  .VI  rliseriminator  has  been 
derived  p.ia'Vion.s]-.  ^  1  ».  19  and  it  has  l>ecn  .shown  that 
>ine  ol  :!s  forms  nsfs  the  sam  e  ronfignrallon  of  clement.s 
a.s  the  pha.se- lock  FM  dlseriminafor  but  that  in  thi.s 
form  it  requires  n  nonrealizatile  loop  filter.  Thi.s 
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optimum  discriminator  can  be  realized,  however,  If  the 
received  data  can  be  stored  on  tape  (as  Is  done  In  some 
telemetry  systems)  and  processed  by  a  special  comput¬ 
er.  This  type  of  demodulation  Is  not  practical  for  most 
situations.  Nevertheless,  the  study  of  its  threshold  per¬ 
formance  provides  a  useful  bound  on  the  threshold  be¬ 
havior  obtainable  with  other  techniques. 

In  this  paper,  the  threshold  is  evaluated  on  the  basis  of 
the  expected  percentage  of  time  the  discriminator  will  be 
In  an  unlocked  state.  When  the  probability  of  being  in 
the  unlocked  state  exceeds  a  certain  prescribed  small 
value  under  optimized  conditions,  the  discriminator  is 
assumed  to  have  reached  its  threshold.  Threshold  is 
evaluated  on  the  basis  of  attainable  operation  with  a  sig¬ 
nal  that  Is  frequency  modulated  by  a  sine  wave  which  has 
a  maximum  modulation  frequency  f  and  has  a  peak 
frequency  deviation  f^j  .  This  modulation  format  is 
generally  useful  both  for  performance  evjduatlon  and  for 
experimental  verification.  Stationary  gaussian  white 
noise  Is  assumed  as  the  interference. 

The  discriminator-transfer  functions  for  both  the  phase- 
lock  ajid  the  frequency -lock  discriminators  are  optim  ¬ 
ized  under  the  constraint  that  they  be  of  second  order  or 
less.  This  constraint  is  imposccl  for  practical  reasons 
concerning  both  tlie  ajialysis  and  the  implementation. 
However,  because  of  the  stability  requirements  of  these 
feedback  circuits,  little  improvement  is  e.xpected  for 
higher-order  transfer  functions. 

Phase- l.ock  FM  IDi.seriminator 


Ucsci^)tion 

A  block  diagram  of  one  form  of  the  phase-lock  KM  dis¬ 
criminator  Is  shown  in  Fig.  1.  The  received  data  arc 
amplified  by  a  band-pass  filter,  which  has  a  b.andwidth 
B,.f  sufficient  to  pass  the  slgnid  components  rel.ntively 
undistorted,  and  (mter  the  multiplier  of  the  phase-lock 
loop.  When  the  discriminator  is  operating  in  its  'Tocked- 
on"  mode,  the  sine  wave  generated  by  the  voltage- 
ccntrollotl  oscillator  (Vt'O)  tracks  the  phase  of  the  in¬ 
coming  signal;  the  multiplier  output  provides  the  phase 
correction  informal  ion  necessaiy  to  keep  the  VCO  track¬ 
ing  properly.  I'he  loop  filter  which  has  a  transfer  func¬ 
tion  F(p)  removes  much  of  (he  noise  and  high  fre<|Uency 
components  from  the  multiiilior  output  and  hence  reduces 
their  deleterious  effect  on  the  tracking  operation. 


Fig.  1— Block  diagram  of  phase-lock  FM  discriminator. 


the  statistics  of  the  amplifier  input  and  output  are  iden¬ 
tical  except  for  the  appropriate  gain  factors,  and,  hence, 
the  SNR  at  the  loop  input  is  the  same  as  that  at  the  AGC 
amplifier  input.  If  a  band-pass  limiter  Is  used  Instead, 
the  input  and  output  SNR  differ  only  slightly  and  are 
related  by  !'<  (SNR)out  =  « (SNR)jn  where  7i/-l  -  ff  s  2. 
If  (SNR)j^  »1  .  we  have  o  v  2,  and  if  (SNR)  «1  , 
we  have  o  ?-  jr/4.  The  noise  statistics  change  consid¬ 
erably  on  passing  through  the  band-pass  limiter,  how¬ 
ever,  and  care  must  be  used  in  applying  this  SNR  re¬ 
lationship. 

The  ideal  AOC  amplifier  is  used  In  this  analysis,  al¬ 
though  the  results  are  .shown  to  apply  closely  to  thebond- 
pass  limiter  as  well.  If  the  output  power  of  the  AGC 
amidifier  Is  denoted  by  I’.p  ,  then  the  output  signal 
amplitude  1'.  is  related  to  the  output  noise  power  by 
e2/2  +  =  p,p  .  Since  Input  and  output  SNR  are 

identical  for  the  amplifier,  the  sign.al  amplitude  is  ex¬ 
pressed  by 


K  =  /~2P.j,/v/rVr 


(1) 


Because  of  the  nature  of  tlie  i>ltase-lock  loop,  the  esti¬ 
mate  of  the  signal  phase  <ilitainc'l  liy  the  feedl)aek  system 
must  be  made  in  real  time.  A  fixed  delay  T  is  usuallv 
tolerable,  how  ever,  in  obtaining  the  final  freciiiencv  e.sti- 
mate  at  the  receiver  output.  Therefore,  post-detection 
filtering  (transfer  function  F,xl  (!') )  is  u.suallv  nece.ssary 
to  obtain  tlie  optimum  frequency  estimate  when  the  T 
sec.  delay  i.s  .allowed  and  the  filter  constraints  .arethere- 
bv  relaxed.  The  operation  of  the  post-dctectlon  filter, 
however,  does  not  affect  the  discriminator  threshold 
and  is  of  only  minor  concern  in  this  imper. 

It  is  common  practice  to  in  c'  ede  the  phase-lock  bx.p 
with  an  .'\GC  amplifier  having  <  onstani  power  output,  or 
a  band-pa.ss  limiter  (whicli  also  has  constant  power  out- 
putl.  These  eiements  restrict  the  dynamic  range  of  the 
loop  input  and  cause  the  loop  gain  to  be  dependent  onlv 
upon  the  input  .SNH.  If  .an  ideal  ACrC  .amplilier  is  u.scd. 


where  r  ^  l/(SNR)j„  . 

I.Inearized  Equivalent  Circuit 

The  fitst  step  in  tin-  analysis  of  the  phase-lock  loop  is 
to  eomimte  the  low-frequeney  terms  that  are  present  In 
the  multiplier  output.  For  this  purpose,  let  us  express 
the  signal  eomponent  at  the  output  of  the  AGC  amplifier 
as 

E  sin  [  a'„l  •  ./.{tf  ] 

where 

.  I't  j  f  (i.  t  'ip 
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the  phase  modulation  produced  by  a  frequency  modulat¬ 
ing  message  waveform  f  (t)  .  Define  the  VCO  output 
as 

s(t)  =  A  cos  +  <^(t)  3 

and  the  noise  at  the  multiplier  input  as 

n(t)  =  N(t)sin[  cc^t  +  V'(t)]  • 

Then,  by  referring  to  Fig.  1  ,  the  multiplier  output 
X  (t)  can  be  expressed  as 

2x(t)/A  =  E  sin  [0(t)  -  ^(t)]  +  N(t)  sin  [(f'(t)  -  0(t)] 

=  E  [sin  €(t)  +  0^(t)  ] 

where  terms  in  the  frequency  range  2f  have  been 
neglected,  the  phase  error  is  denoted  by  e  4  0  -  0, 
and 

0^^(t)  ^  [N(t)/E]sin  [0{t)  -  0(t)  ]  . 

The  first  term  in  this  expression  represents  a  phase 
correction  term,  mid  the  second  represents  interference. 

When  the  discriminator  is  ojierating  above  threshold, 
the  phase  of  the  VCO  is  closely  tracking  the  phase  of  the 
signal,  1.  e.  ,  i  (1)  ~  0  (t)  •  In  this  paper,  the  signal  and 
noise  components  are  considered  to  be  statistically  in¬ 
dependent.  Thus,  when  the  discriminator  is  operating 
above  threshold,  the  power  spectral  density  of  0^(i) 
can  be  e.xpressed  approximately  as  the  low-frequency 
portion  of  the  convolved  VCO  and  noise  siiectra;  i.  c.  , 

i:V.  (f)  '  C.-(f)  *  G^(f)  , 

^  n 


where  C.^  (f)  ,  O  •  (0  .  ^'n***'  ‘■‘’pnesent  the 

power  spectra  of  0„(t),  s(t),  and  n(l)  ,  re.spectively . 


The  phase  correction  term  K  sin  c  l  an  be  aiiproxiniatod 
by  the  piecewise  linear  curve  shown  in  Fig.  2.  In  the 
region  of  phase  error,  I  £  |  1  the  multiplier  output  is 

approximated  by 

2x(t)/A^  k[£(I)  ♦  !■;[<.',  (tl  -<;■  (t)  ]  Gl) 

where  0  i,*>  ‘  0  and  can  tic  termed  llu  e(|uival- 

enl  phase*  input  I  hi's  la.st  c'prcssion  li’.ids  to  the  lin- 
cari/i’d  ci,ai\ ali-nt  circuit  for  the  phase- lock  loopshown 
in  Ki.:  'k 

This  circuit  is  equivalent  to  the  actual  r.onlineart  ircui: 
III  that  thi'  loop-filter  input  is  identical  with  that  givmi 
in  i.’l,  and  hence  the  outputs  of  the  f'.'.o  circuits  arc  the 
same.  Notice  that  in  tliis  context,  the  I'-rni  0ntt)  ''an 
lie  considered  the  e(|uivaleiil  input  phase  noise,  it  can  be 
shown  to  lie  a|ipro,xiniatclv  the  same  a-;  th<  phase  noise 
term  that  nouhi  ai'pea  r  at  the  output  of  a  liandpass  iim- 
iter  at  high  input  .''NP  f>in).  Hence,  the  us-  of  a  iiand- 


PHASE  CORRECTION 


Fig.  2— Phase  correction  term  for  the  phase -lock  loop. 


eouivaleni  phase 

INPUT 


Fig.  .1  Linearized  equivalent  circuit  for  the  pha.sc  -lix’k 
loop. 


pass  limiter  at  high-input  SNK  is  not  expected  to  per¬ 
form  differently  than  with  the  AtK.'  amplifier.  At  low 
SNH.  the  output  Si'K  of  the  l)and[ia.ss  limiter  is  ir/-l 
times  the  input  SNH,  or  'inly  1.0  db  bebiw  the  Input  SNK. 
Hence,  the  results  eomp'ited  on  the  ba.'^.i.s  ol  an  AGti 
amplifier  .should  eorres|ioiid  closely  to  those  for  a  tiand- 
pass  limiter. 

The  transfer  function  of  the  linearized  elrcuit  (equival¬ 
ent  phase  input  to-freqneney  estimate)  is 

7fp)  ^ 

♦p  (P)  1  '  gp  F  (p)/p 

where  7  ip)  and  4’  (I'l  are  the  i  ,i)ilaee  transforms  of 
fft)  and  0p(t)  .  respectively,  and  we  define 
g,^,  -  gAI'./2,  the  equivalent  gain.  Hence,  the  phase- 
transfer  funelion  (etiiiivalenl  phase  input  t'l-phase  e.sti- 
mate  ,  transfer  fiineli"n  and  the  transfer  fiindion  of 
equivalent  fieirie'vv-to-frei)ueney  estimalel  is  express¬ 
ed  liv 


H  (I'l 


•h  'pi 
4,,  (1)1 


g,.  F  (p)  p 
‘  g,.  Ffjii'p 


•ifiere  0  ipi  1^  the  io'iiiei  transform 


ft  I 


(2.) 
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Two  types  of  phase-error  components  are  present  in  the 
loop— transient  error  and  noise  error.  The  transient 
phase  error  is  defined  as  tlie  error  in  (p  which  is  de¬ 
pendent  only  upon  the  input  signal;  it  has  the  transform 
Et(P) 

Ej(p)  =  $(P)[1  -  H(p)]  (4) 

The  phase-noise  error  is  defined  as  the  error  in  the 
phase  estimate  which  is  dependent  upon  the  input  noise, 
and  it  has  the  transform  E_  (p) 


Threshold  is  evaluated  on  the  basis  of  a  signal  which 
has  a  sinusoidal  phase  modulation  of  the  form 

<t>{t)  =  (V^a^  “-a^ 

where  f^  £  f the  maximum  modulation  of  frequency 
and  fjj  is  the  peak  frequency  deviation.  The  transient 
phase  error  with  this  modulation  is  a  sinusoid  of  fre¬ 
quency  f^  having  a  peak  value 

S  =  V^a)  I  1  <®) 


En(P)  =  4-jj(P)H(p)  (5) 

Both  (4)  and  (.‘5),  of  course,  are  based  upon  the  assump¬ 
tion  that  the  linearized  analysis  is  a  good  approxima¬ 
tion.  These  expressions  are  used  in  the  next  section 
to  evaluate  the  discriminator  performance. 

The  dynamic  (frequency-tracking)  range  of  the  dis¬ 
criminator  Af  is  a  fundamental  parameter  of  the  loop 
and  is  defined  as  the  maximum  steady -state  frequency 
shift  that  can  be  produced  by  the  VCO.  The  maximum 
average  multiplier  output  occurs  with  e  =  v/2  and 
has  the  value  =  EA/2.  Notice  that  on  the  basis 

of  the  linear  model .  this  maximum  frequency  shift 
occurs  for  e  =  1 .  The  corresponding  maximum  VCO 
frequency  shift  is  Af  =  EAgF(0)/2  =  ggFfO). 

The  product  F  (0)  represents  the  dc  loop  gain,  and 
it  is  convenient  to  define  the  quantity  g^  =  ggFfO) 
for  use  in  the  next  section. 

Threshold  Effect 

In  the  piecewise  linear  approximation  to  the  phase  cor¬ 
rection  term  of  the  actual  phase-lock  loop,  the  correc¬ 
tion  is  linear  with  respect  to  phase  error  for  |  e  |  <  1 . 

If  (  e  I  should  exceed  1,  the  linear  beh.avior  no  longer 
holds— the  phase  estimate  becomes  highly  distorted, 
and  the  loop  can  be  unstable.  In  this  paper,  the  dis¬ 
criminator  is  considered  to  reach  its  threshold  when 
the  input  SNU  decreases  to  the  point  where  the  proba¬ 
bility  of  having  ]  e  |  >1,  computed  on  the  basis  of  the 
linearized  analysi.c,  exceeds  a  prescribed  small  proba¬ 
bility.  More  specifically,  threshold  is  defined  to  occur 
when  the  rms  value  of  the  total  phase  error  oj  =  1/2. 
A  certain  degree  of  arbitrariness  is  inherent  in  any 
definition  of  threshold. .and  this  one  i.s  no  exception.  It 
is  difficult  to  specify  at  just  what  probability  of  losing 
lock,  or  how  far  below  the  linear  behavior  of  the  output 
SNR  curve,  threshold  should  be  .said  to  occur.  However, 
some  measure  of  the  arbitrariness  is  removed  in  a 
later  p.araprai)h  where  the  probaiiililv  of  losing  lock  is 
computed  as  a  function  of  the  input  .SNR.  The  curve  of 
probability  of  lo.=  ing  lock  vs.  input  .SNR  which  isobtained 
gives  insigiit  into  the  variation  of  the  discriminator 
performance  as  the  input  SNR  vanes  slightly  above  or 
below  its  threshold  value.  In  types  of  FM  discrimina¬ 
tors  wliicii  track  the  phase  or  frequency  of  the  input 
signal  the  percent  time  out  of  lock  is  perhaps  a  more 
meaningful  measure  of  perform.ance  than  the  output  SNR 
because  near  threshold  large  errors  in  the  estim.ate 
occur  in  bursts.  In  the  time  interval.s  between  bursts, 
the  output  noise  level  is  relatively  small,  and  the  accu¬ 
racy  in  the  message  estimate  is  relatively  good. 


The  value  of  peak  transient  error  used  for  the  thres¬ 
hold  determination  is  the  maximum  value  of 


e.  for  f  £  f 
t  am 


The  phase  noise  has  a  mean  square  value 


(f)  I  H(ja)^)  l^df 


-OO 

Thus  the  total  rms  phase  error  is 


,1/2 


(V) 


can  be  computed  from  (6)  and  (7);  this  quantity  deter¬ 
mines  the  threshold. 


The  closed-loop  transfer  function  considered  in  this 
paper  is  constrained  to  have  the  form 


HtP) 


_ 1  +  =tP 

1  +  b^p  +  bgp^ 


(8) 


Appendix  A  shows  that  the  approximately  optimum 
form  of  this  transfer  function  for  large  deviation  ratio 
is 


1  +  ap  +  a  p 

This  transfer  function  is  approximately  optimum  for  the 
deviation  ratios  of  main  interest  in  this  paper 


The  loop  filter  required  for  the  realization  of  this 
transfer  function  is  the  RC  filter  depicted  in  Fig.  4. 
The  closed-loop  attenuation  characteristic  |H  ()(>>)  |  is 
plotted  in  Fig.  •'>.  The  pe.ak  transient  error  occurs  for 

f  =  f  with  this  filter, 
a  r.i 

It  l.s  shown  in  Appendix  A  that  the  input  .SNR  is  related 
to  the  peak  transient  error  and  the  phase  noise  by  the 
expression 


(SNR). 


2  1/2 
’’  r . 

n  t 


rf 


(10) 
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TRANSFER  FUNCTION 


F(p)  = 


1+  ac 

l  +  0*fljP 


Rj=  ARBITRARY  R, »  (opj,-!)  Rj 

C  =  o/R,  9g=  LOOP  GAIN  =  pg  g, 

Oflo  »•' 

Fig.  4— RC  loop  filter  for  the  phase-lock  loop. 


Notice  that  for  large  deviation  ratios,  D  »■  1  ,  the 
threshold  input  SNR  is  proportional  to  _  Equa¬ 

tion  (12)  is  discussed  further  in  the  final  section  of  the 
paper. 

The  probability  that  the  phase  error  will  exceed  the 
threshold  error  |  e  |  =  1  based  on  the  linearized 
analysis,  can  be  computed  as  a  function  of  the  input 
SNR.  If  the  phase  noise  in  the  estimate  is  assumed  to 
be  approximately  gaussian,  then  the  total  phase  error 
is  the  sum  of  a  sinusoid  plus  a  gaussian  term;  the  sta¬ 
tistics  of  this  sum  are  available  from  the  literature^®. 
At  threshold,  the  probability  of  exceeding  |  e  |  =  1  is 
about  0.  043  .  The  probability  oflosing  the  locked-on 
state  is  plotted  In  Fig.  6  under  the  assumption  of  a 
fixed  linearized  closed-loop  transfer  fimctlon  H(p). 


Fig.  5-Closed-loop  attenuation  characteristics  for  the 
phase-lock  loop. 


1/2  1/2 

The  value  of  a  used  is  a  -  c.  /w  D  .  The 

I  m 

quantity  Brf  denotes  the  equivalent  noise  bandwidth  of 
tho  RF  amplifier. 


If  this  expression  is  minimized  with  respect  to  for 
a  prescribed  rms  total  error.it  becomes 


(SNR) 


in 


1.57 


jtD 


1/2 


f 

m 


1/2 


for  =  \/2/5  (r.j,  ,  tr^  =  -v^Ts  ff.j.  ■  Thus  the 
threshold  SNR  is  given  by 


1/2, 

irn  f 


ISNR).^^  =  8.9 


(in 


Consider  that  the  RF  filter  is  single-tuned  and  has  a 
2  db  tiandwidth  2  (D  ^  H  ^m  >  ll*®  equivalent  noise 
liandwidth  is  R,.f  *  t(D  I'^m  ■  The  threshold  ex- 
;iression  (lli  tlicn  becomes 

(SNR)j^.j.  =  8.9  +  1) 


loop. 


r.ccausc  of  nonlinear  effects,  the  probability  of  losing 
lock  is  assumed  to  be  essentially  unity  for  input  .SNR 
l>elow  the  thresliolri  value.  When  the  discriminator 
losc.s  lock,  it  must  go  through  a  rather  complex  lock-ln 
operation  l>eforc  linc.ar  operation  is  resumed.^ 

OulpulSNR 

When  tho  post  -detection  filter  is  proi’.erly  optimized  at 
high-input  SNR  for  a  low-pass  mfKiulalion  f  (t)  ,  the 
combined  transfer  function  of  the  phase-lock  loop  and 
tho  post-rietnetion  filter  has  a  constant  attenuation 
characterisi  i<j  aiifl  a  linear  phase  charaetcristic  for 
'll  '  111!  •  Tilt;  filler  greatly  attenuates  frequency 
components  al>ove  the  modulation  pass -band  if|  >  f^  . 
When  the  probability  of  losing  lock  i.‘-  small,  the  input 
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and  output  SNR  of  the  complete  discriminator  can  be 
shown  to  be  related  by  the  expression 

g 

<SN^)out  =  I  (SNR)j„  ,  (13) 

m 

A  sinusoidal  frequency  modulation  having  a  peak  devia¬ 
tion  =  DL,  has  been  assumed  for  this  calculation, 
d  m  ’ 

Experimental  Results 

A  phase-look  FM  discriminator  has  been  constructed 
with  filter  parameters  based  on  the  theory  presented 
herein:  it  has  the  following  modulation  characteristics; 


Message  bandwidth 

fm 

=  2kc 

Peak  frequency  deviation 

'd 

=  lO.Skc 

Deviation  ratio 

D 

=  5.  25 

Equivalent  noise  bandwidth 

rf 

=  (ir/2)60kc  =  94kc 

This  experimental  discriminator  utilized  a  band-pass 
limiter  rather  than  an  AGC  amplifier.  However,  it  is 
not  expected  to  alter  the  discriminator  performance 
significantly  from  that  produced  by  the  AGC  amplifier , 
because  the  m:un  region  of  interest  is  for  input  SNR  a:  1 
A  measured  curve  of  output  vs.  input  SNR  for  the  dis¬ 
criminator  is  plotted  in  Fig.  7,  A  theoretical  curve  is 


A  measured  curve  of  the  fractional  time  out  of  lock*  is 
shown  in  Fig.  8  along  with  a  theoretical  curve  based  on 
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Fig.  8— Measured  fractional  time  out  of  lock  vs.  input 
SNR  for  the  phase-lock  discriminator. 

Fig.  6.  (Note  that  the  frnctional  out  of-Iock  scale  is 
logarithmic.)  The  agreement  between  these  curves  Is 
quite  good.  The  measured  value  of  the  fractional  time 
out  of  lock  is  0.043  db  for  an  input  SNR  (threshold)  of 
2.  2  db,  which  is  only  0.7  db  above  the  predicted  thresh¬ 
old.  The  measured  fractional  time  out  of  lock  curve 
ch.anges  at  a  rale  greater  than  that  given  by  the  steeper 
linearized  analysis  under  the  assumption  of  gnussban 
noise,  but  does  not,  of  course,  have  a  sudden  disconti¬ 
nuity  at  the  threshold.  A  partial  explanation  for  the 
steeper  slope  can  be  found  by  noting  that  these  measure- 
ment.s  were  made  using  a  bandpass  limiter  and,  ns  a 
result,  the  phase-noise  stalislies  deviate  from  the  g.auss- 
ian  statistics  that  have  been  as.sutned.  If  the  limiter 
effect  is  t.aken  into  account  in  computing  the  curve  in 
Fig.  B,  the  theorcliral  curve  would  be  considerably 
steeper  in  the  region  above  threshold  and  agree  more 
closely  with  the  mca.sured  data  of  k'ig.  8. 


Fig.  T  Mi’.isured  curve  of  output  vs.  input  SNR  for  the 
pha.se -lock  discriminator. 

plotted  in  the  same  figure  which  obeys  (13)  for  iniwt 
SNU  above  the  thre.shold  value  and  drops  to  zero  at  and 
below  threshold.  The  threshold  input  SNR  computed 
from  (11)  is  1.1  db.  There  is  good  agreement  between 
the  aieasurcd  and  the  theoretical  curve.  It  can  be  seen, 
however,  that  the  output  SNR  decreases  noticeably  below 
the  linear  curve  :\t  injnit  SNR,  2  db  above  the  threshold 
v.aluc;  at  threshold,  the  output  SNR  decreases  to  +  21  db. 


*The  (r.actionnl  time  out  of  lock  is  determined  by  mea¬ 
suring  the  fractional  time  in  which  Large  bursts  of 
noise  appear  in  the  loop  output.  These  noise  bursts 
are  at  least  an  order  of  magnitude  larger  than  the 
noi.se  level  between  bursts  and  hence,  are  readily 
scp.arated  from  this  ambient  ntiise. 
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Frequency-Lock  Discriminator 


Description 


The  frequency-lock  FM  discriminator  Is  shown  In  block 
diagram  form  In  Fig.  9.  As  the  figures  shows,  the  In- 

INPUT 


Fig.  9-BIock  diagram  of  the  frequency-lock  FM  dis¬ 
criminator. 


coming  signal  plus  additive  gaussian  white  noise  is  first 
amplified  by  a  band-pass  amplifier  of  bandwidth  suf¬ 
ficient  to  pass  the  signal  with  low  distortion,  and  is  then 
mixed  with  5  si.iusoid  generated  by  the  VCO.  When  the 
discriminator  is  trael' '  •.i  ptoperly,  the  instantaneous 
frequency  :jf  il'i  ’‘fT- r'U-  closely  tracks  the  instanta¬ 
neous  frcquC'ir,  of  Iho  sign.nl  with  a  fixed  IF  frequency 
off-set,  fjp.  Ii  •he  dov!  '  :  :tio  is  greater  than  unity, 

the  condition  o'  "ro:.t-  .  'i  rest  in  this  paper,  then  the 
residual  frequency  deviation  of  the  signal  after  mixing 
can  be  considerably  le.'i  ■  than  th"?  of  the  RF  trans¬ 
mission  Because  of  irdecacn  in  the  deviation,  an 
IF  amplifier  bandwidth  less  than  the  RF  bandwidth  can 
be  used,  and  much  of  the  noise  can  be  removed  before 
the  signal  reaches  the  FM  discriminator  inside  the  loop. 

The  ideal  FM  discriminator  used  in  this  loop  generates 
a  voltage  output  proportional  to  the  instantaneous  fre¬ 
quency  of  its  input.  This  discriminator  provides  a  cor¬ 
rection  voltage  which  is  dependent  upon  the  frequency 
tracking  error.  The  loop  filter  F(p)  serves  to  eliminate 
much  of  the  frequency  noise  from  the  discriminator 
while  it  p.asses  the  modul.ation  comi»onen^s  with  low  dis¬ 
tortion;  tlie  filter  output  is  the  estimate  f(t)  of  the  modu¬ 
lation.  A  post-detection  filler  fpj(P)  is  used  to  obtain 

the  optimum  estimate  of  the  modulation  with  a  time  delay 
T,  i.e.  ,  f  (t  -  T).  As  was  true  with  the  phase-lock 

loop,  the  post -detect  ion  filter  does  not  affect  the  dis¬ 
criminator  threshold. 


Linearized  Equivalent  of  the  Frequency-Lock  Discrim¬ 
inator 

The  analysis  of  the  frequency-lock  discriminator  begins 
with  the  computation  of  the  components  in  the  mixer  out¬ 
put.  It  is  convenient  for  this  purpose,  to  express  the 
signal  and  VCO  output  waveforms  In  the  forms: 

s(t)  =  E  sin  [tu^t  +  <f<t)J 


and 


s(t)  =  A  sin  [(£0^  +  +  ^(t)] 


where  ^ 

0(1)  =  J  f(p)dp 

o 

and 

fV)dp. 

“O 

If  we  again  represent  the  noise  as 

n(t)  =  N(t)  sin  [w^t  +  if'(t)] . 

then  by  referring  to  Fig.  9,  the  mixer  output  x.„(t) 
can  be  expressed  by 

2xjp(l)/AE  =  cos  [<Ujpt  +  (^(t)  -  0(t)] 

+  cos  [U'jpt  +  ll'(l)  - 

=  cos  [uijpl  -t  ^(t)  -  0(t)]  +  n^ft) 

where 

Ogll)  ^  [N(1)/K]cos  [wjj,l  +  il>(l)  -  0(t)]  . 

Terms  in  the  frequency  range  f  +  f^p  are  neglected. 
The  first  term  In  this  expression's  the  frequency  cor¬ 
rection  term;  the  .second  is  an  interference  term  caused 
by  the  noise  Input.  The  residual  frequency  deviation 
in  the  correction  term  can  be  defined  as 

1,(1)  =  f(t)  -  f(l)  . 

Consider  that  the  noise  n(t)  is  gaussian;  then  the  noise 
term  in  the  mixer  outinit  is  also  gaussian.  When  the 
discriminator  is  operating  above  threshold  the  fre¬ 
quency  modulation  and  the  estimate  of  this  modulatton 
are  approximately  equal .  Hence  the  spectrum  of 
n2(t)  is  approximately  equal  to  the  low  frequency  por¬ 
tion  of  the  convolution  flg(f)*Gn(0 

The  mixer  output  is  pa.ssed  through  an  IF  filter  having 
a  symmetrical  (high  Q)  single-tuned  response  ami  a 
.1-db  bandwidth  ft,,  .  The  Ir.ansfer  function  of  this 
filter  is 


LOCKHEED  AIRCRAFT  CORPORATION 


MISSILES  ond  SPACE  DIVISION 


2  43  rU-2 


'  1  +  (p  -  j0Jjj,)/7rBjp 


^  1  +  (p  +  jWjp)/irBjp 

This  transfer  function  can  also  be  expressed  in  the 
forms 

FjpGcu)^ 


=  F^(jw  -  jGjjp)  +  F^(ja)  +  jcojp)  (14) 

The  function  F^  (j  oj)  is  the  low-pass  equivalent  of  the 
IF  filter.  ^ 

The  signal  component  at  the  IF  output  is  distorted  by 
an  amount  which  is  dejiendent  upon  a  complex  relation¬ 
ship  between  peak  value  of  the  residual  frequency  devi¬ 
ation  fj.p,  the  maximum  modulation  frequency  fm  • 
and  the  IF  bandwidth  Bjp  .  BaghdadylS,  however, 
has  shown  that  a  quasi-stationary  response  is  obtained 
with  a  single-tuned  high  Q  fitter  if  4fjjjfj.p  «  In 

the  quasi-stationai y  approximation,  a  sinusoid  of  in¬ 
stantaneous  angular  frequency  ojj(tl  on  passing  through 
a  filter  having  a  transfer  function  iFipnoj))  exp 
is  amplitude  distorted  by  iFjp  [jc»)j(t)]  f  and  phase 
distorted  by  0[iWj(t)]  .  In  this  treatment  of  the  fre¬ 
quency-lock  discriminator  the  IF  bandwidth  is  consid¬ 
ered  to  be  large  enough  so  that  quasi-stationary  analysis 
is  a  good  approximation,  and  the  signal  component  in  the 
IF  output  is  repre.Hetiled  by 


iAE|Fjp[jw^(t)]|co8 


t  + 


e'0‘>’,p)]dii| 

(1.0) 


where  a)^(t)  ^  2)tf^_(t)  is  the  residual  angular  modula¬ 
tion,  and  O'Ooijp)  is  Ihc  group  delay  of  the  filler, 
0'(ja’)‘^  dO(jui)/diu  .  Therefore,  the  main  effect  on  the 
residual  frequency  modulation  is  a  delay  of  O'Oeiii.’) 


•sec. 


Denote  Ihe  IF  output  noise  component  as 

If  the  residual  phase  modulation  of  the  sign.al  is  denoted 
by  0p(t) ,  then  the  frequency  noise  com|)oncnt  In  the 
discriminator  output  is  the  derivative  <if  the  phase  noise 
■  Hy  •I’O  use  of  a  phasor  diagram  this  phase  nouse 
can  be  evaluated  as 

rf-n  larr'  [N,p(t)  K].sin  [i'j^.(t)  -  rf,^(t)] 

'  IpC)  --  r''"*) 


then  the  Fourier  transform  of  0„(l)  ,  denoted  as 
is  approximately  equal  to 

^nO^^’)  -  .A(^(jce)  FaOw)  (16a) 

where  (j£u)  is  the  Fourier  transform  of 

7^„(t)  =  [N(t)/E]  cos  [!f(t)  -  0(t)]  . 

This  last  expression  shows  the  phase  noise  to  be  tlie 
same  as  a  hctcrorlyned  version  of  the  mixer  output 
noise,  namely,  Wa(t)  .  which  has  been  filtered  by  a 
low-pass  equivalent,  F;|(ja))  ,  of  the  IF  filter.  How¬ 
ever,  the  assumption  that  this  statement  is  based  upon, 
namely,  f^(t)  '  0  ,  is  not  completely  accurate. 

Nevertheless,  it  can  be  seen  that  when  Ihc  discrimina¬ 
tor  is  locked  on,  the  signal  component  in  the  mixer 
output,  cos  r(.L)jj,.t  +  i;!,.(i)j  .  is  of  small  enough  band¬ 
width  to  pass  the  IF  relatively  undistorted.  Hence,  the 
additional  phase  modulation  i/>j.(t)  ,  which  appears  in 
(16),  is  not  large  enough  to  change  greatly  the  (lower 
spectrum  of  (/>r,(t)  from  that  which  would  be  obtained 
from  (16a)  . 


Under  the  conditions  slated  above,  ihtreforf, 
put  of  the  ideal  discriminator  is 


'r  -  ''’'<*‘‘’1F>1  ^  2^-df- 


where  k  is  Ihc  gain  of  the  discriminator. 


the  out- 


The  apiiro.ximate  Fourier  iranslorm  of  this  term  is 

kj[.7(.i‘'’)  .7(i‘<.')]  il  Nj^(.k')j  F,^(ju')  (Kia) 

when'  7(.iw)  ,  ,7(iw)  are  Ihe  Fourier  transforms  of 
f(l)  and  l’(l)  ,  r('.s|)ecl  ively .  I'he  a|)|iro.\lmation  is 
made  that  the  dominant  effect  ot  F,  (jtu)  on  Ihe  resid¬ 
ual  frequency  deviation  term  Is  to  c.ause  the  delay 
h'fju’jp)  as  in  (ir,) 

Thus  by  reference  to  Fig.  !)  and  (lOa)  it  c'an  lie  .seen  that 
Ihe  Fourier  Iran.sform  ol  the  frequency  esl i mate  is 
approximatelv  equal  to 

O'^Oa.')  sF  (iw)  F(ja')  [  7(iuj)  -  7(j,.i’)  -  7^^  (ju;)] 

where  7|^|,(iu.') jCV.fju,')  is  tiu'  transform  of  fn,.(l)  , 
the  <  rjnivaU'iiI  fn'cjui'nrv  noisf  input..  VmU.T  rooflilions 
of  lorkcMi-on  opi  ralion,  Ihi'  Infjurnrv  rstimnto  ohtaimMi 
b\  di.scri ?ni nator  (  an  b('  s(’(  n  to  bo 

ih('  snmo  :if  ihni  provirirci  tiv  the  lim  ari/rH  f'rjui vnl(*nt 
('ircMiit  shown  in  Ki;:,  10.  Tlv  -Insi  ft- luop  iransb  r 
function  of  this  f-ircnu  is  ''Hfi*  )  '*7(,1''')  7(V’-)  wherf' 


n{\.) 


A 

7  (Uci 


I  ^{|o  )K(  lo.) 

“g^(V)F(i'..) 


if  |N,j,.(t)l  K 


ami  the  cf’uivnicnf  u.iin 


ku 


A!  SNR  above  threshold,  the  assumption  that 
|N|p(f)i  '  ^  *■'  0  valid  one.  It  ean  be  seen  that  if  lh< 

residual  modulation.  dt .  is  smali. 


As  was  trur  uifb  tie  lin  k  looji.  two  f>l 

errfir  coinpoiii  ct.'-  ai*  |)r'‘d'e  ''fl  in  tlie  > 
proresi;  —  t ra H'' ifiit  f-rror  nifl  rvu.-^i  i  rroi-,  Tlif- 
transient  «rror,  uhich  n 'xiilts  fion;  a 
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Fig.  10-Llnearized  equivalent  circuit  of  the  frequency- 
lock  loop. 


frequency  modulation  having  the  transform  .!7(jtu)  . 
has  the  transform  . 

=  .9(j“)'[l  -  H(ju))]  (17a) 

The  output  frequency  noise  has  the  transform 

-  .7  H(jw)  (17b) 


Threshold 


In  an  "ideal"  FM  discriminator,*  the  outjHil  SNH  begins 
to  decrease,  substantially  more  rapidly  than  the  first 
(tower  of  the  intuit  SNR  when  the  intwt  SNR  goes  as  low 
as  *  10  db.  In  the  frequency-lock  discriminator,  the 
threshold  effect  occurs  when  the  SNR  at  the  input  to  the 
ideal  discriminator  decreases  below  10  dh,  and  signal 
su|)|)ression  effects  begin  to  be  of  signific.ancc.  This 
IF  SNR  can  decrease  to  its  threshold  value  from  two 
causes;  (1)  The  noise  (lowcr  at  the  in|nit  to  the  fre¬ 
quency-lock  discriminator  ran  increase  relative  to  the 
signal  [lower.  (2)  The  residu.al  frequency  deviation  in 
the  mixer  output  can  increase  to  the  (loint  where  the 
signal  com(ionent  falls  outside  the  IF  (lass-band  and  is 
attenuated  relative  to  the  noise  and  otherwise  distorted. 
The  objective  of  the  o()tiinizalion  that  is  (terformed  in 
this  (laiier  is  to  adjust  the  IF  filter  bandwidth,  the  l«o|i 
gain,  the  loop-filter  transfer  function  so  that  these  two 
causes  of  threshold  occur  at  as  low  an  in[Kit  SNR  ns 
(possible. 

The  o|itimi7.ation  is  to  be  (lerformed  tiv  niiniini/.ing  tae 
ms  total  frequency  error  Ojj  by  proiierly  choosing 
the  loop  transfer  function  coefficients.  The  IF-filter 
3-dti  bandw  idth  Bjp  is  then  given  the  smalle.st  (mssible 
vaiiic  to  (lerinil  significant  distortion  of  the  sign.i.1  com- 
(loncnt  (IF  frequency  correction  com(ionent)  to  begin  .at 
tiiicshold.  The  signal  com|ioncnt  in  the  IF  is  considered 


to  require  a  mininuitn  bandwidth  of  approximately 
2(ffjj  +  /2  o-fT)  for  low  distortion  when  the  maximum 
modulation  frequency  is  fm  •  Notice  that  for  sinusoi¬ 
dal  frequency  error,  the  peak  error  deviation  is  /2 
times  the  rms  frequency  error.  Threshold  is  said  to 
occur  when  the  IF  bandwidth  has  Lhi.s  value,  and  the  IF 
output  SNR  decreases  to  10  db. 

The  frequency  modulation  of  the  signal  is  assumed  to 
have  the  form  f^j  cos  ui^t  where  uijj  s  01^^  .  When  the 
discriminator  is  locked-on  and  the  linearized  analysis 
(17a)  is  accurate  the  i.eak  transient  error  f,,,  has  the 
value 


^tp  =  ^dl 

The  rms  frequency  noise  error  can  be  found  from  (17b) 
to  be 


^  f  G,  (f)  |H(jcu)|2df  (19) 

n  -"o  ‘ne 

where  Gf^^(f)  is  the  power  spectral  density  of  f„e(t)  . 
the  equivalent  frequency  noise  input. 

The  closed-loop  linearized  eijuivalent  transfer  function 
of  the  frequency-lock  discriminator  is  constrained  to 
have  the  form 

H(p)  --  g/(1  +  bj  p  +  bg  (1^) 

A()pendix  B  shows  that  the  approximately  optimum  form 
of  this  transfer  function  is 

H(|))  - - S - - - 

1  +  1)1)  +  bp  /to 

‘  '  ill 

The  loop  filter  used  to  obtain  this  transfer  function  is 
shown  in  Fig.  il.  The  [leak  transient  frequency  error 
(18).  now  occiii'.s  when  the  modulation  frequency  is  at 


Rj  r  arbitrary 

R.  =  R2(-fo-') 

C  =  0i  /  R^ 

L  =  (  02  /  oq  )  Ra 


LOOP  GAIN  go  =  «z  »  I'**?  "m 
Oo»  l/W  B|F 
0,  :  ( !♦  go)  b 
02=  (  I  go)  b/ta)m 


*  Broad -bandw  idth  (niisc  counting  or  Fo.stcr-Secl\  tviies 

of  FM  discriminators  arc  good  atiproximations  to  an  Fig.  il  RI.C  lorq,  filter  for  the  fre.piency-loek  dis- 

idcal  FM  discriminator  crimiontoi 
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Its  niaxlmum  f  .  With  this  transfer  function,  the  peak 
transient  and  noise  frequency  errors  are  related  to  the 
Input  SNR  by  the  equation 


(SNR)  = 


B 

rf  f^  tp 


This  relationship  is  obtained  from  (19)  by  setting 

b  =  f.  /2^'\  . 

tp  tp  m 

the  value  obtained  by  solving  (18)  for  b  with  f^  =  f^ 
Now,  tor  a  given  (SNR).^,  the  rms  total  error 


Is  minimized  with  respect  to  by  setting 


where  B  ,  is  the  noise  bandwidth  of  the  RF  filter.  Bv 
rf 

substituting  in  (24)  the  e.xpresslon  for  (23),  we 

obtain 


lOirf  I  .  27rf 

(snr)^=^^iM2D)  / 


If  we  define  (SNR)^B^^/27rf^^ x,  then  this  equation  can 

be  rewritten  as  D  -  1)^^ x  ,  a  relationship  plot¬ 

ted  in  Fig.  12.  As  can  be  seen  from  Fig.  12,  if 

-too  c - 1 - ^ . . 1 - 1 


ASmPTOTE  D«t 
1*5 


—  ASYMPTOTE  D>*l 

■  1 

■  L  1  1  1  III 

DEVIATION  RATIO 


By  using  (21a)  the  rms  total  frequency  error  then  can 
be  expressed  by  (see  Appendix  B) 


Fig.  12-Normalized  threshold  input  (SNR)t  as  a  func¬ 
tion  of  the  deviation  ratio  for  the  frequency-lock 
disc  rlmlnator. 


’"fT  “  m  °  B^f(SNR)j^ 


The  minimum  IF  bandwidth  consistent  with  this  fre¬ 
quency  error  is  then 

B,p  =  2(f^W2o„) 


2f^  V  -  (2») 


[B^f(SNR).rJ 


where  the  threshold  values  (as  yet  unknown)  of  n^j.  and 
Input  SNR  are  used  in  (23). 

The  threshold  value  of  input  SNR  occurs  when  the  IF 
SNR  reaches  10  dh.  Since  the  equivalent  noise  band¬ 
width  of  the  single-tuned  IF  filter  is  rBjj-/2,  the  thresh¬ 
old  Input  SNR  (SNR).j.  is  expressed  by 


D  >>  1,  then  x  ~  (.lOD)  '  is  a  good  approximation.  For 
D  «  1,  the  normalized  SNR  x  approaches  a  minimum 
value  of  5.  This  curve  allows  the  delermination  of  the 
threshold  SNR  for  .an  arl)ltr.arv  RF  liandwldth  B 

rf 

A  relation  of  even  more  interest  is  the  curve  of  thresh¬ 
old  SNR  as  a  function  of  the  deviation  ratio  0  when 

B  -  =  7r(D  +  1)1  ,  its  minimum  value  for  low  distortion, 
rf  '  m 

Then  the  threshold  can  be  obtained  from  Fig.  12  by  the 
relation  (SNR).j,  -  2x/(D  +  1),  and  the  result  Is  shown  in 

Fig.  13.  As  can  be  seen  from  the  figure,  the  threshold 


joml  fm  oiscrjwinatw  |_ 

-  FREQUENCY -LOCK  OISCRMlNATON 


_ 1 


ASYMPTOTE  0-;l 
(SNR),  -OKO 


(SNR)jp  =  (SNR).^  rf  -  10 


DEVIATION  RATIO  0 


(SNH).^  -  dttBjj. 


Fig.  13  Throsholfi  (.SNillq-  a.s  .n  function  of  deviation 
ratio  for  th<'  fri  rnicnc v-lock  di.scrimlnator. 
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is  approximately  10  db  for  D  s  1,  and  it  decreases  ap- 
-3/1 

proximately  as  9 . 6D  for  D  >  1. 

Experimental  Measurements 

An  experimental  frequency-lock  FM  discriminator  has 
been  constructed  in  accordance  with  the  theory  pre¬ 
sented  in  this  paper,  and  its  performance  has  been 
measured.  The  discriminator  was  designed  to  demod¬ 
ulate  a  wide  deviation  FM  signal  having  the  following 
parameters: 

Maximum  modulation  frequency  f^  =  l.SKc 
Peak  frequency  deviation  f ^  -  1  Me 

Deviation  Ratio  D  =  667 

A  single-tuned  RF  filter  is  used  which  has  a  3-db  band¬ 
width  of  2  MC;  hence,  =  :r(l  Me)  =  ^-(0  +  l)fj^- 

Flgure  14  shows  the  measured  fractional  time  out  of 
lock  for  the  discriminator  as  a  function  of  the  input 
SNR.  The  theoretical  threshold  value  of  input  SNR  based 
from  Fig.  13  is  -7.0  db.  This  prediction  agrees  well 
with  the  measurements  which  gave  a  5  percent  time  out 
of  lock  at  input  SNR  between  -6.. I  and  -7.0  db. 


Fig.  14  Mea.sured  value  of  the  fr-actional  time  out  of 
lock  vs.  input  SNR  for  the  freqiieney-loek  di.scrim- 
inator. 


The  outpul  SNRforIhe  disc riniinalor « a.*;  al,«o  measureil 
as  :i  function  of  the  input  SNR.  Because  of  extraneous 
circuit  noise,  hum,  eU-.  .  the  highest  output  SNR  mea¬ 
sured  was  al  db.  The  output  SNR  did  not  decrea.se  .sub¬ 
stantially  below  this  value  until  the  input  SNR  was  re¬ 
duced  to  a  threshold  v.alue  of  between  -5, .4  and  -6.5  dt., 
iust  as  expected  from  the  other  measurement.  The 
threshold  effect  was  very  prominent.  When  the  input 
SNR  was  decreased  to -S  ",dl'.  the  output  SNR  decreased 
sharply  to  approximately  unitv 


Maximum  Likelihood  FM  Discriminator 
Description 

As  mentioned  in  the  introduction,  the  maximum-likeli¬ 
hood  FM  discriminator  can  take  the  same  form  as  the 
phase-lock  loop.  However,  in  that  form,  it  requires 
a  nonrealizable  filter.  As  a  particular  example,  let  us 
consider  that  the  frequency  modulation  has  a  low-pass 
power  spectral  density  Gf(f)  of  the  form 

Gf(f)  =  fj/4f  for  I  f  I  s  f 
f'  '  d'  m  '  '  m 

(25) 

=  0  I  f  I  >  I 

'  m 

1/2 

where  f(j/2  “  is  the  rms  frequency  deviation,  and 

fjj^  is  the  maximum  modulation  frequency.  The  inter¬ 
fering  noise  n  (t)  is  considered  to  be  band-limited 
gaussian  white  noise  having  a  constant  spectral  density 
Gf,(f)  "  N(j  in  the  signal  pass  band.  Under  these  con¬ 
ditions,  the  maximum-likelihood  k’M  discriminator  can 
be  .shown  to  have  the  form  shown  in  Fig.  15.  The  loop 
filler  F„(p)  is  nonrealizable  because  it  corresponds 
to  a  filter  impulse  response  which  is  nonzero  for  t  ^  0. 


RECEIVED  DATA  FREQUENCY  ESTIMATE 

E  tin  (wot*  ^  (t)) 

♦n(t)  a?(t) 


Fig.  15  Ma.xlmum-llkellhood  FM  discriminator 


An  approximate  com|)Ulcr  realization*  of  this  estimator 
is  show  n  in  Fig.  16  This  l  ompiiler  o|)'^  r  ite  s  on  recorded 

*The  maximum  likelihood  estimate  can  also  be  obtained 
byusingan  a  (xisteriori  probabilitv  computer.  Theeom 
puter  of  Fig.  16  i."-'  discussed  here  because  it  most  re¬ 
sembles  the  form  oftlie  phase-lock  loop.  Roth  realiza¬ 
tions  required  the  use  of  a  (on  plcx  tv|ie  of  waveform 
generator 


LOCKHEED  AIRCRAFT  CORPORATION 


MISSILES  and  SPACE  DIVISION 


RECORDED  TIME  SEGMENT 
or  RECEIVED  DATA 


E  Iln  (Wot  ♦  ^(0)  r,()w)-rFo{Jwl«'^“'^ 

♦  n(l)  * 

f,{t-T) 


Fig.  16-An  approximate  computer  realization  of  the 
maximum-likelihood  FM  discriminator. 


time  segments  of  the  received  data  t  sec.  long.  It 
uses  a  complex  type  of  waveform  generator  which  gen¬ 
erates  a  large  number  of  modulating  function  which  are 
samples  t  sec.  long  of  a  waveform  having  a  power 
s|)ectruin  Gf(f)  as  in  (25).  Obviously,  the  number  of 
waveforms  to  be  generated  is  intolerably  large  for  any 
but  the  moat  simple  classes  of  modulation  functions. 


The  linear  filter  used  in  Fig.  16  is  an  arbitrarily  good 
approximation  to  the  optimum  F^(p)  and  is  realizable 
because  it  allows  an  arbitrary  delay  of  T  sec.  in  the 
estimation  process.  The  computer  operates  by  com¬ 
paring  the  waveforms  f((t)  at  the  filter  output  with  a  de¬ 
layed  version  of  the  test  modulating  waveform  fj(t).  The 
pair  that  matches  Is  the  maximum  likelihood  estimate 
and  that  fj(t)  is  the  Identical  waveform  (except  for  a 
delay)  as  is  obtained  by  the  estimator  of  Fig.  15. 


To  compute  the  threshold  of  the  maximum  likelihood 
estimator,  we  make  the  identical  set  of  calculations  as 
were  made  for  the  phase-lock  loop,  except  that  the 
closed-loop  transfer  function  is  now  of  the  form 


4 

♦  (Jw) 


II  (Jw) 


i-^r.f(f)/rN^ 

1  ^  r  (:,(f)/f^N 
s  f  o 


where  I’g  =  F)''/2  ,  the  signal  power,  and  Gf(f)/f“  1“ 
the  po^ver  spectral  density  of  the  phase  mixiulation 
if«(t).  Define 


X  -  2f  /n  ,f.  (SNR). 

m  rf  d  in 

where  (SNR)jj^  =  Pg/2R|.fN^  and  is  the  equiv¬ 

alent  noise  bandwidth.  Then  the  closed-loop  transfer 
function  simplifies  to 


-  1/(1  .  \f“)  If  I  -  f^ 

('  i  f  !  >  f 

m 


(26) 


The  parameters  of  this  filter  change  with  the  (SNR)(n  . 
However,  in  this  paper  where  threshold  is  the  most 
important  effect,  the  filter  Is  fixed  and  optimized  for 
the  threshold  value  of  (SNR)j_  .  Thus,  the  filter  time 
constant  Is  chosen  X  =  2f^;Bj.f(SNR)T., 

As  in  the  previous  computations ,  threshold  Is  computed 
on  the  basis  of  a  sinusoidal  frequency  modulation, 
f  (t)  =  fjj  sin  tUjj,  t .  The  peak  transient  phase  error  Is 
Cj  =  2(X/D)/(1  i  \/D  ).  The  mean  square  phase 
noise  error  Is 


<jI  =  r  G  (f)  I  H  (jcu)  I  ^df 

-oo 


where  (f)  is  the  equivalent  in])ut  phase  noise  spec- 

trum  and  Ga  =  1/2  B  (SNR),  in  the  low-frequency 
1 1  in 

region  where  H(ju>)  >  0  The  mean  square  phase  noise 
can  be  evaluated  as 


2 

(j 

p. 


dl 


(1  ' 


2 

xn 


=  .  Vi__  ^ 

1  +  X/D- 


Define  the  threshold  lo  occur  when  the  rms  total  error, 
(Trj<  reaches  1/2.  The  erjuation  defining  threshold  is 

ft^/2  I  -r^  =•  (1/2)^ 


and  can  be  written 


1  _a 

4  i  +  a 


1/2 


Ian 


h 

(2  0)“ 


4  2 

where  a  -  X/D  This  expression  can  he  solved  for 
1)  us  a  function  of  a  and  then  the  threshold  (SNR)|n 
can  be  obtained  by  relating  it  to  a  .  The  result  is 
shown  In  Fig.  17.  For  deviation  ratios  D  >  1  ,  Jigood 
appro.ximation  is  (,SNIl)-|  ^  1.  27/D 


Fig.  17  Threshold  (.SNR)jn  .is  a  function  t,f  the  devia¬ 
tion  ratio  for  the  m.nximum-HkcHhood  estimator 
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Comparison  of  the  Discriminators 

Each  of  the  three  types  of  discriminators  has  been 
evaluated  on  the  basis  of  Its  threshold  value  of  Input 
SNR  for  a  prescribed  type  of  frequency  modulation. 

When  optimum  post-detection  filtering  Is  used  and  the 
discriminators  are  operating  at  SNR  well  above  their 
thresholds ,  the  relationship  between  the  output  and  in¬ 
put  SNR  is  the  same  for  each  type.  Thus,  it  is  quite 
meaningful  to  compare  the  performance  of  the  discrim¬ 
inators  on  the  basis  of  their  thresholds. 

hi  Fig.  17,  the  threshold  input  SNR  Is  plotted  as  a  func¬ 
tion  of  the  deviation  ratio  for  each  of  the  discriminator 
types.  In  each  example,  the  equivalent  noise  bandwidth 
Bj-f  Is  assumed  to  be  7r(D  +  l)fm>  the  minimum 
value  consistent  with  low  distortion  of  the  signal 
waveform.  The  +  10  db  threshold  SNR  of  the  conven¬ 
tional  or  "ideal"  discriminator  is  shown  in  the  figure 
for  comparison. 

At  small  deviation  ratios  D  <  1  the  thresholds  of  the 
two-phase  coherent  dlscrlminator.s  (phase-lock  and 
maximum-likelihood)  actually  decrease  with  decreases 
in  the  deviation  ratio.  The  reason  for  this  effect  is  that 
the  maximum  phase  excursion  caused  by  the  modulation 
is  less  than  1  radian  in  this  region.  Thus ,  transient 
errors  are  of  little  significance  in  determining  the 
threshold,  and  the  noise  can  be  sharply  filtered.  Of 
course,  for  small  deviation  ratios  and  low  input  SNR 
the  output  SNR  is  so  small  as  to  be  of  limited  usefulness. 
Furthermore  the  deviation  caused  by  drifts  in  the  oscil¬ 
lator  center  frequencies  may  have  to  bo  taken  into 
account  with  this  type  of  operation. 

The  frequency-lock  discriminator  shows  no  improve¬ 
ment  over  the  ideal  KM  discriminator  for  small  devia¬ 
tion  ratios.  The  explanation  for  this  behavior,  of  course, 
is  that  the  IF  filter  has  a  minimum  bandwidth  of 
Bif  =  2ffii  ,  the  same  as  the  RF  bandwidth  B^f  . 
Hence,  the  IF  filter  docs  not  remove  any  significant 
amount  of  noise  from  the  signal. 

At  deviation  ratios  11  >  1  all  three  types  of  discrim¬ 
inators  provide  a  reduced  threshold  input  (SNR).  As 
expected,  the  maximum  likelihood  FM  discriminator 
has  the  lowest  threshold  SNR  for  all  values  of  the  devia¬ 
tion  ratio;  at  high  deviation  ratios,  it  h.as 
(SNR)t  -  1.27/D  .  The  reason  for  this  superiority'  is 
that  the  discriminator  is  not  conslraiiied  to  real  time 
operation,  and  the  loop  filtering  can  thereby  be  much 
more  effective.  However,  as  has  been  .slated,  the  com¬ 
plexity  of  the  waveform  generator  required  of  this  type 
of  discriminator  m.ikes  it  unsuitable  for  all  but  the  most 
elementary  classes  of  modulation  functions. 

H  is  shown  in  Fig.  18  that  the  ph.asc-lock  discriminator 
is  siqierior  in  perform.ance  to  the  frequency -lock  dis- 
crimin.alor  for  small  and  moderate  values  of  the  devia¬ 
tion  ratio  (D  <  25).  This  behavior  can  be  made  plausi¬ 
ble  by  jw)inting  out  that  the  phase-lock  disciimin.ator  is 
a  phase-coherent  device,  while  the  frequency -lock  dis¬ 
criminator  tracks  instantaneous  frequency. 


0.1  10  10  100  1000 
DEVIATION  RATIO  0 


Fig.  IS-Comparison  of  the  threshold  input  SNR  for  the 
discriminators  as  a  function  of  D. 


In  other  words,  when  tlie  discriminator  is  locked-on 
the  VCO  output  is  phase  coherent  with  the  signal  and 
is  noncoherent  with  the  noise.  As  a  result,  the  noise 
power  in  the  multiplier  output  has  decreased  relative 
to  the  signal  component  by  a  factor  of  1/2  from  its 
value  before  the  mulliidieation  operation.  (This  change 
in  noise  level  is  analogous  to  that  found  in  the  synchro¬ 
nous  detection  of  AM  signals.)  Therefore,  the  phase 
coherence  of  the  phase-lock  loop  can  be  an  Important 
advanl.tge  and,  evidently,  is  of  dominating  Importance 
for  deviation  ratio.s  D  <  25.  At  larger  values  of  the 
deviation  ratio,  the  frequency -lock  discrimin.ator  equals 
or  surpasses  performance  of  the  phase-lock  loop. 

This  effect  can  be  attributed,  at  least  In  part,  to  a 
basic  difference  in  the  threshold  errors  of  the  two  dis¬ 
criminators.  The  threshold  phase  error  for  the  phase- 
lock  loop  is  fixed,  I  f  I  =  1  is  used  in  the  linearized 
an,'ilysis,  whereas  threshold  frequency  error  for  the 
frequency -lock  discriminator  Is  dependent  upon  an 
adjustable  quantity,  the  IF  b.andwldth.  The  IF  bandwidth 
can,  of  course,  be  optimized  for  a  given  type  of  mod¬ 
ulation,  and  In  this  sense  provides  a  degree  of  freedom 
tliat  does  not  exist  in  the  phase-lock  loop.  It  Is  inter¬ 
esting  to  note  from  (2.3)  that  If  the  IF  bandwidth  is  fixed, 
perhaps  at  its  optimum  value  for  D  =  1.  then  the  fre¬ 
quency  errors  liecomc  of  dominant  Importance  as  D 
Is  increased.  In  this  circumstmee,  the  threshold  in¬ 
put  SNR  i.s  proportional  to  for  large  1)  ,  the 

s.ame  relationship  as  obtained  with  the  phase-lock  loop. 
By  optimizing  the  IK  bandwidth  a.s  is  done  in  this  paper, 
the  threshold  input  SNR  can  vary  in  proixjrtlon  to 
and.  as  a  result,  allows  the  frequency-lock  discrimin¬ 
ator  to  surpas.s  the  performance  of  the  phase-lock  dis¬ 
criminator  at  large  deviation  ratios. 
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The  results  of  this  paper,  however,  do  not  permit  the 
conclusion  that  the  frequency-lock  technique  Is  superior 
to  the  phase-lock  technique  at  large  deviation  ratios 
where  higher-order  transfer  functions  are  allowed. 

The  analysis  becomes  substantially  more  complex, 
however,  as  the  order  of  the  transfer  function  is  In¬ 
creased  to  third,  or  higher  order.  Furthermore,  con¬ 
straints  must  then  be  placed  on  the  closed-loop  trans¬ 
fer  function  in  order  that  the  loop  filter  required  be 
realizable.  The  absence  of  poles  In  the  right  half¬ 
plane  of  the  closed-loop  response  Is  Insufficient  for 
rellzabllity  of  the  loop  filter. 
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Appendix  A 

of  the  Phase-Lock  l,oop  Transfer  Function 

In  this  appendix,  the  phase  lock  loop,  closed-loop  trans¬ 
fer  function  Is  assumed  to  have  the  form. 


M(p)  .  --.-1.' 

1  ^  hjp  f  bgp 

The  parameters  of  this  function  are  to  be  optimized  so 
as  to  minimize  the  rms  total  phase  error  in  the  modu¬ 
lation  estimate.  Irfiin  (6)  the  transient  phase  error  can 
be  ev.aluatcd  as 


‘  t  f 


Hfjw 


a) 


2  .  ■! 
‘^a  ^2  > 


1,2  4 
*’2  “a 

2,  2 
a;  b, 
a  1 


1/2 


The  t,  .nsient  error  is  largest  for  f  ^  f  the  max 
imum  modulation  freipiency,  and  this  value  will  be 
assumed  henceforth.  For  large  deviation  ratios  and 
I  £  I  r  1  ,  this  expression  ran  li.-  simplified  to 
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These  conditions  of  large  deviation  ratio  and  |  ej  |  be¬ 
low  the  threshold  value  are  of  major  importance  in  this 
paper,  and  it  is  for  these  conditlons'^that  the  closed-loop 
transfer  function  is  optimized. 


The  mean  square  phase  noise  error  can  be  computed 
from  (7)  and  found  to  be 


^l-L 


G.  (f)|Haw)rdf  =  ■ 

-oo 


o,<0) 


1  + 


(28) 


where  G  ^  (f)  G  .  (0)  for  deviation  ratios  D  >  1.  It  can 

K  K 

be  shown  from  (2)  that  the  power-spectral  density 


G^^<«)  =  Kf(SNR)i„ 


and  the  mean  square  noise  error  (28)  is 


o2=[2aB^,(Sm),„] 


-1 


Hence  the  Input  SNR  Is  related  to  the  transient  and 
noise  phase  errors  by  the  equation 


Trf 


(SNR)  = 


Thls  result  Is  the  same  as  (10)  In  the  text. 


Appendix  B 


where  B  ,  is  the  equivalent  noise  bandwidth, 
rl 

Thus,  the  total  mean  square  phase  error  is  of  the  form 

2  2,-^2 

Ot  =  €/2  + 

and  has  the  representation 


"  2B^f(SNR)j„ 


^2*  a 


"m‘’l'’2 


Define  the  quantities 

Then  this  expression  can  be  rewritten 


o 


2 

T 


«) 


^1^2 


It  is  easily  seen  that  the  optimum  value  of  a  is  bounded 
by  0  <  o  <  .  A  good  approximation  to  the  approxi¬ 

mate  minimum  of  is  obtained  with  a  =  and  this 
value  of  n  will  be  used.  The  optimum  value  of  is 
then  /3j,  and  the  resulting  close-loop  transfer  function 
is 


Optimization  of  the  Freauencv-Lock  Transfer  Function 

The  closed-loop  transfer  function  of  the  frequency-lock 
discriminator  is  considered  to  have  the  form 


H(p)  = 


1 

1  +  bjP  + 


When  sinusoidal  frequency  modulation  at  frequency  f^^ 

is  applied  to  the  signal,  the  peak  transient  frequency 
error  can  be  found  from  (18)  to  be 

ftp=fd''-»«“a)l 


=  f 


1.2  2  ^  .2  4 
b,u)  +  b-ci> 
la  2  a 

""T  2\2  ,  ^2  2 

-^“a) 


"1“'8 


The  maximum  transient  error  occurs  for  fa  *  the 
maximum  modulation  frequency,  and  this  modulation 
frequency  will  be  used.  For  large  ratios  of  the  fre¬ 
quency  deviation  to  the  transient  frequency  error,  the 
situation  of  greatest  Interest,  this  expression  reduces 
to 


(f  /f.)-bfen^  +  hlJ^ 

tp  d'  1  m  2  m 


The  mean  square  value  of  the  frequency  noise  error  can 
be  computed  from  (19)  and  has  the  value 


o 


2 

f 


n 


on 


Gj  (f)lH6u))|^df 
ne 


H(p)  = 


1  +  ap 


2  2 

1  +  ap  +  a  p 


The  peak  transient  error  (27)  now  has  the  value 


r,  2  2 

f  =  D  a  0)  I 

t  m 


where  G^  (f)  is  the  power  spectrum  of  the  equivalent 
ne 

frequency  noise  input. 


Of  (f)  =  f^/ 
no 


2Brf(SNR),^ 
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Thus  G.  (f)  Is  a  parabolic  power  spectrum.  Evaluation 
*ne 

of  this  Integral  gives  the  result 


2 

167r='B^j(SNR)j^bjb2 


This  expression  is  symmetrical  with  respect  to  /Jj^  and 
and  Is  minimized  by  setting  P,  =  l.e.,  b„  =  b./ai  . 

A  1  a  1  in 

Hence  the  optimum  form  of  H(p)  Is 


H(p)  = 


]  +  h, 


1 _ 

2  , 

•  op  ,u 
1  m 


If  we  now  define 


^  =  Vm*  ^2  =  ^>2“™-  «  ="4/  [2Brf<SNR),„]  • 


then  the  mean  square  total  frequency  error, 


is 


„2  _fd  /.2  .2\ 

\  -  2  ^2) 


The  resulting  mean  square  error  Is  then 

This  error  Is  minimized  If 


and  thus  the  minimum  mean  square  frequency  error  Is 


27rf 


B^,(SNR)in 


1/2 


The  square  ruot  of  this  quantity,  gives  (22)  of  the 
text. 
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